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We report the in-plane and out-of-plane heat conductivity of the antiferromagnetic spin-orbit 
induced Mott insulator Sr 2 lr 04 with J e g = 1/2. Our data reveal clear-cut evidence for magnetic 
heat transport within the Ir 02 planes which provides the unique possibility to analyze the thermal 
occupation and scattering of J e g = 1/2 pseudospin excitations. The analysis of the magnetic 
heat conductivity yields a low-temperature (T < 75 K) magnetic mean free path Z mag « 32 nm, 
consistent with boundary scattering. Upon heating towards room temperature, the mean free path 
strongly decreases by one order of magnitude due to thermally activated scattering of the pseudospin 
excitations. The latter reveals that the coupling of these excitations to the lattice is radically different 
from that of S = 1/2-excitations in cuprate analogs. 

PACS numbers: 71.70.Ej, 44.10.+i, 66.70.-f 


The physics of iridium oxide materials has recently 
moved into the focus as these materials realize a plethora 
of novel quantum magnetic phases based on, e.g., the 
square, honey-comb, and hyperkagome lattice types [1—5] 
with J e ff = 1/2 pseudospins. This includes the sought- 
after possible realization of a quantum spin-liquid with 
peculiar elementary excitations [6, 7]. The magnetic heat 
conductivity is considered an important tool to probe 
quantum spin and topological excitations [6] , as it is sen¬ 
sitive to both the thermal occupation and the scattering 
of such quasiparticles. In the past years, this sensitiv¬ 
ity has been exploited extensively for probing the ele¬ 
mentary spin excitations in many different 5 = 1/2 low¬ 
dimensional quantum magnets [8-17]. For materials with 
a strong spin-orbit coupling (SOC) and J e g = 1/2 pseu¬ 
dospins, such as the iridates, it remains however unclear 
whether heat transport can be used to probe the pseu¬ 
dospin excitations at all, because the strong SOC is likely 
to cause strong scattering due to phonons. Accordingly, 
successful experiments on iridate materials are lacking, 
apart from pioneering attempts [18]. 

One of the up to the present best studied iridate ma¬ 
terials is the compound Sr 2 lr 04 which is a spin-orbit in¬ 
duced Mott insulator [2] with localized electrons on the 

Ir 4 + 

ions in a J e g = 1/2 state. The material possesses a 
very similar structure as La 2 Cu 04 , i.e. a square lattice 
of Ir 4+ ions is formed by corner-sharing IrC >2 plaquettes, 
where adjacent Ir02-planes are separated by SrO layers 

[19] . A strong antiferromagnetic exchange of the order 
J ~ 0.06 eV couples the J e g = 1/2 pseudospins giv¬ 
ing rise to two-dimensional (2D) magnetic excitations as 
is revealed by resonant inelastic x-ray scattering (RIXS) 

[20] . Sr 2 lrC >4 orders long-range antiferromagnetically at 
Tn ~ 224 K, where a weak ferromagnetic moment occurs 
due to canting of the IrOg octahedra [3, 21, 22]. 


In this Letter, we report the in-plane and out-of-plane 
heat conductivity of Sr 2 Ir 04 . We observe a highly un¬ 
usual in-plane heat conductivity with anomalous tem¬ 
perature dependence at T < T)v which is absent for the 
out-of-plane direction. This is incompatible with phonon 
heat transport and evidences magnetic heat transport 
within the Ir02 planes. Thus, our data reveal the first 
example for magnetic heat transport in a J e g = 1/2 
compound. We analyze the low temperature magnetic 
heat conductivity K mag in terms of a Boltzmann-type ap¬ 
proach and extract the magnetic mean free path Z mag . At 
low temperature, l mag is limited by boundary scattering. 
Upon heating to room temperature, l mag decreases by 
an order of magnitude, indicating temperature-activated 
scattering. This scattering process is the dominating one 
at elevated temperature and can be assigned to strong 
magneto-elastic coupling. 

The growth and characterization of single crystals 
of Sr 2 lrC >4 have been described in Ref. 1. The crys¬ 
tal dimension in our experiment was 0.51 x 0.87 x 
0.17 mm 3 where the shortest edge is the crystallographic 
c-direction. The thermal conductivity k has been mea¬ 
sured with a home-made device in a four-point configu¬ 
ration using a chip resistor as heater and a thermocouple 
to measure the temperature gradient parallel to the ab- 
planes [23]. Due to the limited size of the thin plate-like 
single crystal it was impossible to perform a four-point 
measurement along the c-direction. Nevertheless, a two- 
point measurement was possible. For our setup we care¬ 
fully investigated the differences between the two- and 
four-point configuration. For the two-point configuration 
we found, that above ~ 150 K the temperature depen¬ 
dence of the heat conductivity is reproduced correctly 
with the caveat of an uncertain absolute value. There¬ 
fore, we consider data from the two-point measurement 
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FIG. 1. Heat conductivity k of Sr 2 lr 04 measured in the ab 
plane {n a b) and along the c direction (k c ). Purely phononic 
fits to the heat conductivity are shown as lines. The fit to 
the ab direction is given by the values an = 2.8 W/Km, 
bi = 530 W/m and To = 18 K. Similarly, the phononic fit 
to the c direction with a c = 1.43 W/Km, b c = 160 W/m. 
For the ab direction extreme fit parameters were used to 
find limits of the phononic temperature behavior marked as 
the shaded area. The upper bound is found by the values 
ao = 4.811 W/Km, bo = 284.4 W/m and To = 0, the lower 
bound by a 2 = 1.2 W/Km, &2 = 700 W/m and To = 25 K. 


only to search for anomalous temperature dependence in 
the vicinity of T/v- 

Fig. 1 shows the measured heat conductivity of Sr 2 lr 04 
along the ab and c directions, n a b and k c , respectively. 
From resistivity measurements [24, 25] and applying the 
Wiedemann-Franz-law it follows that Sr 2 lrC >4 has a low 
contribution of electrons to the heat conductivity which 
is three orders of magnitude lower than the measured k. 
Therefore, we neglect the contribution of the electrons 
and consider the total heat conductivity in ab direction 
as the sum of a conventional phononic and a potential 
magnetic contribution, whereas k c is purely phononic. 

The in-plane heat conductivity K a b exhibits a peak at 
low temperature (~ 12 K), which is followed by a steep 
decrease that slightly levels off at around 75 K. At fur¬ 
ther increased temperature a broad step around T/v is 
observed and the curve almost saturates close to room 
temperature. This temperature dependence is incompat¬ 
ible with canonical phononic heat conduction. In a sim¬ 
ple approach the heat conductivity is proportional to the 
specific heat cy, the velocity v, and the mean free path l 


of the heat carriers: 

k ~ cyvl. (1) 

In the case of phonons as heat carriers, the velocity and 
the mean free path are approximately constant at low 
temperature and k follows the temperature dependence 
of the specific heat. At higher temperatures, umklapp 
scattering becomes important which reduces the mean 
free path and thus leads to the observed low-temperature 
peak. This process depends on the number of excited 
phonons and leads to l oc 1/T. Thus, at high tem¬ 
peratures where cy approaches the Dulong-Petit con¬ 
stant, the phononic heat conductivity is approximated 
by K phon oc 1/T [26]. 

The leveling off at ~ 75 K and the broad step-like fea¬ 
ture near T/v are clearly inconsistent with this expected 
temperature dependence. Two completely different sce¬ 
narios are conceivable for explaining this unexpected be¬ 
havior. It is possible that enhanced scattering of phonons 
occurs due to critical magnetic fluctuations near Tat. In 
fact, a dip structure near T/v is often found in antifer¬ 
romagnetic materials like in MnO [27] or in C 0 F 2 [28]. 
However, such critical fluctuations are unlikely to affect 
Kphon near 75 K, i.e. far away from T/v. Moreover, the 
phonon scattering due to magnetic fluctuations typically 
affects the heat conductivity isotropically even in layered 
systems [29]. Therefore, we carefully inspected the heat 
conductivity parallel to the c-axis in the vicinity of T/v- 
At temperatures higher than ~200 K, the k c curve is 
absolutely featureless and fully described by the afore¬ 
mentioned 1/T-law as is indicated in the figure. Thus, 
we can clearly rule out a phononic origin of the anoma¬ 
lous behavior in n a b- On the other hand, the anomalous 
behavior can also arise from a 2D magnetic heat conduc¬ 
tivity within the Ir02-layers which adds to Kphon, i.e., n a b 
results from the sum of phononic and magnetic contribu¬ 
tions while k c is purely phononic. Indeed, such magnetic 
heat transport is frequently observed in low-dimensional 
5 = 1/2 quantum magnets [8-17]. Thus, we conclude 
that in Sr 2 lr 04 a magnetic contribution to the heat con¬ 
ductivity is present in n a b- 

Having established this main experimental finding, 
namely the first observation of magnetic heat conduc¬ 
tivity in a Jeff = 1/2 system, we move on to its quantita¬ 
tive analysis by extracting the magnetic mean free path. 
A phenomenological T _1 approach is used to model the 
phononic heat conductivity at high temperatures (cf. the 
solid line in Fig. 1). To estimate the uncertainty of 
the phononic background in n a b , we performed extreme 
phononic fits for determining lower and upper bounds as 
is indicated by the shaded area [30]. 

In the temperature range between T/v and ~ 50 K the 
measured k exceeds the expected phononic heat conduc¬ 
tivity remarkably, corroborating our conclusion of a sig¬ 
nificant K mag . We subtracted the phononic fit from the 
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FIG. 2. Magnetic heat conductivity K mag of Sr 2 lrC >4 resulting 
from the difference between k and the phononic fit (dots). 
The line is the low-temperature fit with constant magnetic 
mean free path l mag and a constant shift w from Eq. (2). The 
uncertainty from the phononic fit is spanned by the shaded 
area. 


measured k and obtain K mag as shown in Fig. 2. K m ag in¬ 
creases with increasing temperature up to ~ 150 K with 
a maximum value of about 1.8 W/Km, and decreases for 
higher temperatures. If one considers the coarse generic 
behavior of heat conductivity given by Eq. (1), the low- 
temperature increase arises from the thermal occupation 
of pseudospin excitations (termed magnons hereafter). 
The peak at 150 K and the high-temperature decrease 
of Kmag cannot be related to the maximum of the spe¬ 
cific heat of 2D spin excitations, as the latter is expected 
around 0.6 J/k B ~ 420 K [31] (with J the exchange con¬ 
stant in a 2D 5 = 1/2 Heisenberg model), i.e. at much 
higher temperatures than considered here. Instead, the 
decrease must be primarily related to a temperature de¬ 
pendence of the mean free path since the average magnon 
velocity is unlikely to change strongly at the rather low 
temperatures (T -C J/k Bl with J = 0.06 eV [20, 32, 33]) 
considered here. 

For investigating our result for K mag further, we follow 
an approach that has previously been used to analyze 
Kmag of the S = 1/2 analog La 2 Cu 04 [15], i.e., we extend 
the simple kinetic description of Eq. (1) by accounting 
for possible momentum dependencies in two dimensions 
[15, 17], i.e. k oc f CkVkludk, with Ck = the 

specific heat (ek and n k are the energy and the Bose 
occupation function of the mode k), Uk the velocity and 
/k the mean free path of a magnon with wave vector k. 

In Sr 2 lr 04 , the dispersion £k is a steeply increasing 
function with band maxima at e > 100 meV = fcs-1160 K 
[20]. Thus, at the low temperatures considered here, pri¬ 
marily modes with small momenta are relevant for the 
heat transport. For simplicity, we therefore assume a 
temperature independent mean free path / mag = Ik- We 


approximate the dispersion ek = Ck = y/A 2 + (Hvok) 2 
[15], where A = 0.83 meV is the experimental anisotropy 
gap revealed by ESR [34] and Vo = (6.25 ± 0.4) • 10 4 m/s 
the small-fc magnon velocity extracted from the RIXS 
single magnon dispersion [20]. Thus, considering two 
magnon branches, we get [15]: 


Om/T 

_ , 2v 0 k B T 2 f 2 l~ “ e x dx 

h mag “ mag a 2 c Q 2 m J X V ■* ( e * - l) 2 ’ 

Xq 

( 2 ) 

with Xo = A/k B T and a = 3.9 A, c = 25.8 A the relevant 
lattice parameters [35, 36], and ©m = hvo^/ir/ak B the 
Debye temperature for magnons. 

We use Eq. 2 to model the low temperature K mag by 
assuming a temperature independent / mag for a certain 
temperature range. This corresponds to the physical pic¬ 
ture of dominating magnon boundary scattering. A cor¬ 
responding fit is shown in Fig. 2 where we account for a 
possible offset w due to an uncertain phonon background 
at low temperature [15]. The fit describes the data well 
up to 75 K and yields a low temperature mean free path 
of /mag = 32 ± 10 nm. This value corresponds to ~ 82 
times of the Ir-Ir distance a. 

We now address the apparent deviation from a 
temperature-independent mean free path that becomes 
evident for higher temperatures T > 75 K. In this regime, 
the afore used simple low-temperature Debye approach 
cannot be employed anymore because it fails to properly 
describe the magnetic specific heat c mag (T) at elevated 
temperatures. We therefore use the theoretical result 
of c mag (T) of the 5 = 1/2 Heisenberg antiferromagnet 
on a square lattice [31, 37] with the exchange coupling 
J = 0.06 eV [20, 32, 33] and the above Vo for calculating 
the temperature dependent l mag (T) based on the kinetic 
expression (1) for a 2D system. At low temperature (c.f. 
Fig. 3), Zmag(r) is roughly constant up to ~ 75 K re¬ 
flecting the low-temperature boundary scattering as re¬ 
vealed by the afore analysis. For increasing temperature, 
/mag (T) decreases strongly where the change amounts up 
to an order of magnitude near room temperature. At 
the highest temperature accessible in this experiment, 
/mag seems to saturate close to 1.2 nm, i.e., one order of 
magnitude above the Ir-Ir distance which constitutes a 
natural minimum value of / mag . 

The strong decrease at elevated temperatures clearly 
signals the onset of a temperature-activated scattering 
process. We assume that both the low-temperature 
boundary scattering and the temperature-activated 
process are independent of each other. Following 
Matthiessen’s rule, the mean free path can then be writ¬ 
ten as 


/ 


-l 

mag 



/exp (TVT)V 1 
V A s T j ’ 


(3) 


where the second term is an empirical formula for 
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FIG. 3. Magnetic mean free path l m s.g(T) of S^IrCU (dots) 
resulting from the phononic fit to k along ab. The shaded 
area marks the range spanned by the lower and upper bound 
for the phononic heat conductivity. At low temperatures 
the constant magnetic mean free path is found to be lo = 
(32 ± 10) nm. The dashed line is the phenomenological 
fit to l m ag (T) by Eq. (3) with T * = 482 K and As = 
1.4 ■ 10' K _1 m _1 . The inset shows the magnetic mean free 
path together with the spin-spin correlation length £ (open 
squares) determined by RIXS [36] on a semilogarithmic scale. 


temperature-activated scattering in magnetic heat trans¬ 
port that has been successfully used in one-dimensional 
systems [10, 12, 38, 39] (with T* the characteristic energy 
scale of the temperature-activated scattering process and 
A s a proportionality factor). As can be seen in the figure, 
this formula fits the experimental / ma g(T') quite well. The 
fit yields T* ~ 480 K [40] which should be considered as a 
very coarse estimate of the energy of the most important 
scattering mode. The value roughly lies in the energy 
range of Ir-O-Ir bond bending modes, which have been 
suggested to strongly couple to the electronic structure 
[32, 41], Thus, the primary cause of the temperature- 
activated scattering may be ascribed to the scattering of 
the magnons off these phonons. 

It is instructive to compare the mean free path with the 
spin-spin correlation length £ measured by resonant x-ray 
diffusive scattering [36], see inset of Fig. 3. This quan¬ 
tity is a conceivable natural upper limit for Z mag . In the 
long-range ordered phase below X)v, the spin-spin corre¬ 
lation length is infinitely large and thus unimportant for 
the magnetic heat transport. For higher temperatures, 
£(T) decreases rapidly with increasing temperature but 
remains still more than an order of magnitude larger than 
^m a g• Therefore, it plays only a minor role in limiting the 
magnetic heat conductivity above Tjy, if any. This sug¬ 
gests that the seeming anomaly in n a i > around Tjv (cf. 
Fig. 1) is mostly unrelated to the onset of magnetic or¬ 
dering but rather connected with the growing importance 
of temperature-activated magnon scattering. Note, how¬ 
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FIG. 4. Magnetic heat conductivity K ma g of Sr 2 lr 04 in com¬ 
parison with that of La 2 Cu 04 (reproduced from Ref. 15). 

ever, that a faint change of slope is discernible in the 
semilogarithmic representation of Z mag shown in the in¬ 
set of Fig. 3. 

The strong magnon-phonon scattering evident in our 
data reveals a qualitative difference of the pseudospin 
heat transport of Sr 2 lr 04 and the spin heat transport of 
the almost isostructural and thus closely related S = 1/2- 
system La 2 CuC >4 [15], as is inferred from a direct compar¬ 
ison, cf. Fig. 4. Both the spin wave velocity of La 2 Cu 04 
[42], and the low-temperature magnon mean free path 
(~ 56 nm) of the sample considered in the figure [15] 
are approximately twice larger than those of our Sr 2 lr 04 
sample. Considering these parameters and Eq. (2), the 
almost identical low-temperature increase of both curves 
at T < 100 K can be understood as the consequence of 
dominating magnon-boundary scattering in both cases. 
However, upon further increasing T, a strong suppres¬ 
sion of K mag of Sr 2 Ir 04 occurs while that of I^CuCH 
continues to increase up to room temperature. Appar¬ 
ently, the magnon-phonon scattering in Sr 2 Ir 04 is dra¬ 
matically stronger than that in La 2 Cu 04 , despite similar 
phonon spectra in both compounds [32, 41, 43]. This un¬ 
ambiguously evidences a peculiar and particularly strong 
nature of the magneto-elastic coupling in Sr 2 lr 04 , aris¬ 
ing from the large SOC and the resulting entanglement 
of spin and orbital degrees of freedom [22]. 

In conclusion, our data provide the first experimen¬ 
tal result for low-dimensional J e s = 1/2 pseudospin heat 
transport in an iridate compound. Our data show that 
the magnetic heat conductivity remains a valuable tool 
to probe the generation and the scattering of magnetic 
excitations also in these systems. At low temperatures 
T < 100 K, the magnetic heat transport K mag is dom¬ 
inated by magnon scattering off static boundaries and 
thus comparable with that of 2D S = 1/2 systems. How¬ 
ever, at higher temperatures unusual strong magnon- 
phonon scattering becomes increasingly important, high¬ 
lighting the peculiar nature of the pseudospin moments 
and excitations. 
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